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O^ , Abstract. - We study the electromagnetic coupling and concomitant heating of a particle in 

("^ . a miniaturized trap close to a solid surface. Two dominant heating mechanisms are identified: 

0^ ' proximity fields generated by thermally excited currents in the absorbing solid and time- 

0^ ' dependent image potentials due to elastic surface distortions (Rayleigh phonons). Estimates 

for the lifetime of the trap ground state are given. Ions are particularly sensitive to electric 
proximity fields: for a silver substrate, we find a lifetime below one second at distances closer 
than some ten /im to the surface. Neutral atoms may approach the surface more closely: if 
they have a magnetic moment, a minimum distance of one /^m is estimated in tight traps, the 
M , heat being transferred via magnetic proximity fields. For spinless atoms, heat is transferred by 

inelastic scattering of virtual photons ofi^ surface phonons. The corresponding lifetime, however, 
is estimated to be extremely long compared to the timescale of typical experiments. 



The last few years have witnessed an increasing interest in tightly confining traps of cold 
particles. These devices allow to envisage a broad spectrum of applications ranging from 
single-mode coherent matter wave manipulation and low-dimensional quantum gases [^ g, g, || 
to quantum logical registers g, g]. Since steep trapping fields exist near surfaces, traps in 
their vicinity enjoy increasing popularity. This raises the question at what timescale the cold 
particles in these "surface assisted traps" will be heated up, and how they are coupled to the 
nearby bulk which is typically at room temperature [H, H, H| ■ The question is of primordial 
importance for the above-mentioned applications since the heat transfer to the trap inevitably 
destroys the coherence of the matter waves jg] . 

In this Letter, we outline simple models that allow to compute the lifetime of the trapped 
particle which is limited due to its coupling to thermal excitations of the nearby solid. The 
interaction with thermal blackbody radiation is certainly a candidate for a mechanism of 
heating and decoherence. Estimates given by Wineland and Dehmelt [Q and others H, ^ 
show, however, that this source is negligible for typical trap configurations. This is mainly 
due to the fact that the trapped particles are most sensitive to the field fluctuations at the 
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Fig. 1. - Schematical setup of the experiment. 



resonant trap oscillation frequency (a few MHz at most) which is rather low compared to 
thermal frequencies which are in the THz range. More importantly, the resonant photon 
wavelengths are at least several meters. This means that the particle is always located in the 
near field of its macroscopic environment where the electromagnetic field fluctuations differ 
from the free-space blackbody field [Q Q . The excitations of the solid that give rise to this 
near-field effect come in two species: fluctuating electric currents related to the dissipation in 
the solid (finite electric conductivity), and elastic waves (Rayleigh phonons) that propagate 
along the surface. Current fluctuations generate electric and magnetic flelds above the surface 
("proximity fields") that couple to the particle's charge, spin or polarizability. Surface waves, 
on the other hand, distort the electrostatic image of the particle in the solid and lead to a 
time-dependent image potential. We find that ions are particularly sensitive to proximity 
fields and estimate a typical lifetime of less than a second for distances smaller than 10 /im 
above a metal surface. Atoms, being neutral particles, are less affected by the presence of the 
"hot" surface: for a nonzero magnetic moment, they survive several minutes even at distances 
of a few micrometers. Finally, spinless atoms are completely decoupled from the surface at 
experimentally relevant time scales. 



Model. - We consider a particle in the ground state |0) of a one-dimensional harmonic trap 
(oscillation frequency wt) that is oriented along the unit vector n (see flgjl|) To simplify the 
calculations, we assume that the trap is located above a flat surface whose distance z from the 
trap center is large compared to the size a — {h/{2Mu;t)y^^ of the ground state wave function 
(M is the particle's mass). The interaction potential for the particle is then of the form 

T/(i)--x.F(r,t) (1) 

where x = xn is the displacement of the particle from the trap center and F(r, t) is a fluctuating 
force field at the trap position r. If this force shows fluctuations at the trap frequency, it may 
resonantly excite the particle to the first excited trap state |1). According to second-order 
perturbation theory (Fermi's Golden Rule), this happens with a rate 

ri^o(r) = a-'|(l|x|0)|'^n,n,5i^'(r,r;-wt) (2) 

where (l|a;|0) is a "dipole matrix element" (equal to a for a harmonic trap), and the cross 
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correlation tensor of the force fluctuations is defined by 

+00 
S'^ir',r;u)^ j dt {E,{v' ,t)F,{v,Q))T e'-\ (3) 

~CXD 

where i,j denote cartesian indices. The average is taken in thermal equiUbrium at temperature 
T. Note that the "heating rate" Fi^o also governs the decay of the coherence between different 
trap levels, as is easily shown by deriving a master equation for the particle's density matrix in 
the usual Born-Markov approximation. Decoherence in ion traps was theoretically investigated 
by several authors ||§, [^j |§l ^ |2| Q , particularly in the wake of a recent experiment by Meekhof 
et al. jlj]. In distinction to the present work, refs.||l^, |3[ focussed mainly on the decay of 
off-diagonal density matrix elements. 

Ion heating. - The simplest case is that of a trapped ion. The force in (p^) is given in 
terms of the electric field as F = qE. To compute the electric field fluctuations, we use the 
fluctuation-dissipation theorem that relates their spectral density to the field's Green tensor 
Gij{v' ,v;u)), i.e., the field created at r' by an oscillating point dipole located at r [ [lO[ : 

S'^(r,r;uj) ^ ^——j-j^liRG^j{r,r;uj). (4) 

(In our units, fee = 1.) The Green tensor contains a free-space term Im G^ ' = cu^ / {6Treoc^)Sij 
that describes the black-body field, and a surface-dependent term due to the reflection at the 
surface. For typical trap frequencies, the electromagnetic wavelength is much larger than z so 
that we may use the quasi-static approximation for the Green tensor and get (in SI units) 



167reo z^ e(^) + 1 



Im GW (r , r ; c.) = ^^^^ Im ^^^^^ (5) 



where Sij is a diagonal tensor with Sxx = Syy — ^, Szz = 1- A reasonable model for the 
dielectric function e{uj) in (pi) is that of a Drude metal ll^. For frequencies below the electronic 
damping rate, the electrostatic reflection coefficient in^^) becomes 

e(uj) - 1 
Im = 2ujgeo (6) 

e(iLj) + 1 

where g is the metal's specific resistance at the trap frequency. Finally, we may also take the 
high-temperature limit of the Bose- Einstein factor in (0) . This gives a surface-induced heating 
rate for the trapped ion that follows a power law as a function of distance z 

Typical values are shown in fig .0 where we see that close to a silver surface at room temperature, 
the ion's lifetime is shorter than 1 s as soon as the trap gets closer than about 10 /im. The 
figure also shows that close to the metal surface, the electric proximity field dominates over 
thermal blackbody radiation (the dotted line). We note that this result may be understood 
easily using the analogy between the trapped ion and an oscillating electric dipole that is 
damped both radiatively and nonradiatively. (This analogy has already been exploited to 
interpret decoherence times in early ion trapping experiments, see ref. ItI.) Recall that at zero 
temperature, an excited oscillator decays to its ground state because of photon emission or 
nonradiative energy transfer into the environment (here, the absorbing metal surface) pS. 
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Fig. 2. - Transition rate 0—^1 ('heating rate') for a trapped ion above a Ag surface, plotted vs. 
the trap distance z. Trap oriented perpendicular to the surface (n = e^) with oscillation frequency 
ajt/27r = 1 MHz. Solid line: coupling to electric proximity fields (ion mass M = 40amu, charge 
g = e, and temperature T = 300 K). Dotted line: coupling to the blackbody field. Dashed line: 
coupling to thermal voltage fluctuations in an ion trap, as discussed in [p, pi (endcap distance 1mm, 
circuit resistance K= lf2). The shaded rectangle indicates the size and inverse lifetimes of typical 
experiments |tl4|, [l8[. 



For distances well below the transition wavelength, the decay is dominated by nonradiative 
transfer from the dipole's near field to the metal, as described by eq.(||). At low frequencies 
or, equivalently, high temperatures, the decay rate 1 ^ is dominated by stimulated emission 
and hence proportional to the temperature, see eq.(M). On the other hand, the principle of 
detailed balance implies that the inverse transition 0^1 that we are presently interested in, 
occurs at the same rate in this regime. The ion is hence heated up because thermal energy 
is transferred from the metal surface to the ion's near field in a nonradiative way. (A similar 
reasoning has been presented by Anglin, Paz, and Zurek p7[ , though they focussed on the 
energy transfer from a fast-moving ion into a solid and the concomitant decoherence.) 

We may also compare the heating rate (|7|) to the calculations of Wineland and Dehmelt 
[0 , Lamoreaux [g , and James § who studied the heating due to thermal voltage fluctuations 
across the endcaps of a Paul trap. These authors' results are recovered (up to a geometrical 
factor of order unity) if we replace in (|7|) the quantity gj z with the electric resistance i? of the 
endcap circuitry, and interpret z as the endcap distance: the quantity RT / z"^ then gives the 
power spectrum of the thermal electric field between the endcaps. The present model shows, 
however, that even in front of a single surface, electric proximity fields leak out of the solid 
whose power spectrum increases even stronger with decreasing distance. 

Magnetic proximity fields. ~ Trapped ions and atoms frequently have a magnetic moment 
fj, and are therefore heated by time-dependent magnetic fields. The force derives from the 
Zeeman interaction 



V{r,t) = -fi-Bir,t). 



(8) 



As mentioned in the introduction, magnetic fields are created by fiuctuating currents in the 
solid. Following the seminal work of Lifshitz |19|, the spectral density of these currents is 
proportional to the imaginary part of the dielectric function e(r; w) . It has recently been proven 
that the introduction of fiuctuating currents also provides a consistent framework to quantize 
the electromagnetic field in absorbing and dispersive dielectrics [pOl pi). We use the following 
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Fig. 3. - Transition rate 0^1 ('heating rate') vs. trap distance z for a trapped particle with mass 
M = 40aniu. The trap is oriented perpendicular to the surface and has frequency ut/2Ti = 100 kHz. 
Solid line: coupling to magnetic proximity fields above a silver substrate (magnetic moment /i/(27r?i) = 
1.4MHz/G = 1 Bohr magneton, spin i, specific resistance g = 1.6 x 10~® f2cm llq]). Dashed line: 
same above a glass substrate {g = lO^'^ rjcm). Dash-dotted line: coupling to surface phonons via the 
van-der-Waals interaction {c3/{2-Kh) = IkHz^m^, Ag substrate with huJo = 225 K, Ms — 108 amu, 
T] — 0.75 |P21). The heating rate due to blackbody magnetic fields (not shown) is of the order of 
lO-^'^s-i.^ 



representation of the current Fourier transform j(r; lo) in terms of spatially uncorrelated boson 
operators f (r; lo) pi| 



j(r;a;) = iu^/2h£^)Jme{r^uJ){{r;u;) 



(9) 



For a nonmagnetic solid filling the half-space z < 0, the Biot-Savart law yields the following 
magnetic field cross correlation tensor on the vacuum side z', z > 



S'iiIl,z',R, 



^o^^feo Im £{^) 



U 



1 - e-'*"/^ 87r(z' -I- z) 



(10) 



where R = (x, y) denotes coordinates parallel to the surface and Uj is a diagonal tensor with 
txx — tyy — I and tzz — 1- For simplicity, we focus on a trap oriented perpendicular to the 
surface. Taking gradients with respect to z and z' and putting z' = z, we get the desired 
spectral density of the (z-component of the) magnetic force. Using again the low-frequency 
limit for the dielectric function, Im £{u}) = l/(eoi^£'), and the high-T limit of the Bose-Einstein 
factor, we end up with the following heating rate 



ri^o = 



^,lT (Sf^'-nl) 



(11) 



hujtMg 1287rz3 

This quantity is represented in fig.^ for a trap close to an Ag substrate, and one observes a 
relatively large heating rate of the order of 10~^ s~^ at a distance of 1 /im. A glass substrate 
gives a much smaller heating rate (dashed line) because its resistance is larger. We conclude 
that one has to avoid either metal surfaces or particles with spin if one wants to store atoms 
coherently over timescales longer than, say, a few minutes. 



Heating of neutral, spinless particles. - We now turn to the heating of a neutral, spinless 
particle. One would expect it to be less sensitive to fluctuating electric fields, although these 
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couple to the atomic polarizability a via the Stark shift V — §E^. We have computed the 
corresponding heating rate along similar lines as for the trapped ion and found extremely small 
values (below 10~^^ s~^ even at distances around 100 nm above a Ag substrate). In fact, atom 
heating is dominated by a different effect: the distortion of the surface by thermal oscillations 
leading to a time-dependent image potential. The corresponding force may be easily computed 
using the following effective interaction (strictly valid for a rarefied solid and in the quasi-static 
limit): we integrate a 1/r^ dipole-dipole interaction between the atom and the half-space filled 
by the solid. For a flat interface, one finds the usual van-der-Waals potential —c^/z^. If the 
interface is corrugated, but with an amplitude small compared to the distance z, one gets a 
first-order correction V'^^^ (r, t) of the form p^] 

y(i)(r,i)=I]5(Q;^)wQWexpiQ-R. (12) 

Q 

In this expression, Q is a two-dimensional 'lateral' wave vector (in the xy-p\ane parallel to the 
non-excited surface), R are the lateral coordinates of the trap center, the UQ{t) are the elastic 
displacements of the surface at wave vector Q, and, finally, the coupling coefficients are given 
by .g(Q; z) — —{3c3Q'^/{2z'^))K2{Qz) with K2 a modified Bessel function. 

For simplicity, we restrict ourselves to Rayleigh waves ||2^ whose displacement amplitudes 
are confined to the vicinity of the surface. Their amplitudes may be written in terms of boson 
operators 



hrjTT 



where the dimensionless parameter 77 specifies the decay of the surface wave into the bulk (it 
depends on the ratio between the bulk and Rayleigh wave sound velocities) , N^ is the number 
of surface atoms per quantization area, Mg is their mass, and o^d is the Debye frequency [ p2[ . 
Using thermal Bose-Einstein statistics for the Rayleigh waves, we can compute the spectrum 
of the surface oscillations at the trap frequency. In the result, the magnitude of the Rayleigh 
wave vector is fixed to Q = i^t/^R where wr is the Rayleigh wave sound velocity. Typically, 
the sound wavelength is much longer than the trap distance, and the coupling constant (/(Q; z) 
may thus be approximated by its asymptotic form for Qz <^ 1. The heating rate then equals 



Tel 72n^r]{l + ■q'^) 

hujtujlMM. 



Ti-o - ^^^^^pTiT -n ■ (14) 



As shown in fig.|[ reasonable parameters give a still very small heating rate (below 10^^ s^^ 
at 100 nm distance). Even on a timescale of hours, the spinless atom is thus decoupled from 
the thermal excitations of the solid. We note that in evanescent field traps, the interaction 



with light scattered off surface/bulk impurities may be a dominant heating mechanism. (^ 



Conclusion. - In the vicinity of a solid surface at room temperature, fluctuating electric 
and magnetic fields couple to trapped particles and induce a finite lifetime of the trap's ground 
state. If one wants a lifetime longer than about one second, ions must not be closer to the 
surface than about some tens of micrometers. This means that miniaturized coherent ion traps 
are difficult to realize at room temperature. Magnetic proximity fields are weaker, and atoms 
with spin live for many seconds in the ground state even at distances of a micron. The best 
candidates for long-time storage are spinless atoms because they are nearly insensitive to stray 

(^)C. Henkel, R. Adams, H. Gauck, D. Schneble, T. Pfau, C. I. Westbrook, and J. Mlynek, in 
preparation. 
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fields. Their ground state lifetime is much longer than hours and is mainly limited by surface 
waves that distort the electrostatic image potential. 

The preceding estimates have been obtained from simple models for the trap geometry, 
neglecting the finite height of the trap potential and assuming a homogeneous substrate. An 
obvious extension would be to allow for layered media. Another point is the inclusion of a 
cross-coupling between fluctuating and static trapping fields. This increases the heating due 
to the Stark potential, e.g., because the atom acquires a static electric dipole moment in the 
trap. Finally, in some traps (magnetic or near-resonant optical) the coupling to non-trapped 
internal states of the particle, e.g., hyperfine or magnetic states, may lead to a relevant loss 
rate. 
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